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We report the experimental results of high field ESR measurements in magnetic fields up to 
about 53 T on the quasi two-dimensional triangular lattice antiferromagnet NiGa2S4. From the 
temperature evolution of the ESR absorption linewidth, we find three distinct temperature regions, 
(i) above 23 K, (ii) between 23 K and 8.5 K and (iii) below 8.5 K. The linewidth is affected by the 
dynamics of i?2-vortices below 23 K and one of conventional spiral resonance modes well explains 
the frequency dependence of the ESR resonance fields much below 8.5 K. Furthermore, we find an 
anomaly of the magnetization curve around one third of the saturation magnetization for i?_Lc, 
corresponding to the softening of another resonance mode. These results suggest an occurrence of 
Z-z vortex-induced topological transition at 8.5 K. 

PACS numbers: 75.30.Kz, 75.30.Wx, 75.50.Dd 



Two-dimensional (2D) geometric frustrated systems 
have attracted considerable interests over the three 
decades because they are expected to have a novel ground 
state such as a spin liquid. Among them, the 2D trian- 
gular lattice antiferromagnet (TLAFM) has been studied 
extensively both theoretically and experimentally. Since 
a spin liquid state was theoretically suggested in the 
5=1/2 Heisenberg TLAFM in 1973 Q, a lot of atten- 
tion has been paid to the TLAFM till today. Recent the- 
ories, however, proposed an ordered ground state with 
a 120° spin structure in the 5=1/2 TLAFM with only 
the nearest-neighbor exchange interactions |2| and a spin 
liquid in the 5=1/2 TLAFM with distant neighbor and 
multiple-spin exchange interactions [3j, |4|. Experimen- 
tally, there are only a few samples which are highly 2D 
and show no long range order (LRO) down to extremely 
low temperatures Q. Therefore, it was difficult to clar- 
ify what kind of state is realized actually as a ground 
state of such a strongly frustrated system, and ideal 
compounds having highly 2D characters and isotropic in- 
teractions have been strongly desired to investigate the 
ground state of them. Recently, a new quasi 2D 5=1 
TLAFM NiGa2S4 has attracted substantial interests as 
a candidate with the spin liquid state [g|. This com- 
pound has slabs consisting of two GaS layers and one 
MS2 layer which stack along the c-axis, and each slab is 
separated by a van der Waals gap, resulting in highly 2D 
character. In NiS2 layer, Ni 2+ (3d 8 , 5=1) ions form an 
exact triangular lattice. Despite the fact that strong AF 
interactions are suggested from a large Weiss tempera- 
ture 6» w ~-80 K, it has no LRO down to 0.08 K|, and 
some noteworthy characteristics were observed [fj. The 
neutron diffraction measurements have clarified that the 
AF correlation with a wavevector q — (0.158, 0.158, 0) 
extends to only about 2.5 nm (about seven lattice spac- 
ing) even at 1.5 K, which corresponds to a short-range 
order (SRO) with about 57° spiral spin structure. The 



magnetic susceptibility at 0.01 T exhibits a weak peak at 
about 8.5 K (T*) and a finite value at the low-T limit, 
indicating the existence of gapless excitations. The mag- 
netic specific heat shows broad peaks at about 10 K and 
80 K, and T 2 -dependence with no field dependence up 
to 7 T below 10 K. The T 2 -dependent specific heat in- 
dicates the existence of the fc-linear excitation similar 
to the Nambu-Goldstone mode associated with a break- 
ing of a continuous symmetry. The T 2 -dependent spe- 
cific heat was also reported in some 2D frustrated spin 
systems. These curious experimental results have stimu- 
lated further theoretical and experimental studies on this 
compound in recent years H IS Eil El El El El ■ In 
some theoretical works, the blinear-biquadratic exchange 
interactions were considered to account for such peculiar 
experimental results @,[i,E3|- Ver Y recent NMR, NQR 
and /xSR measurements evidenced the occurrence of a 
novel state where the spins have a MHz scale dynamics 
through a weak freezing below about 10 K Q ■ 

In this Letter, we report the results of ESR measure- 
ments in high magnetic fields up to about 53 T on single 
crystals of NiGa2S4 to clarify the spin dynamics in more 
detail. We have found that the dynamics of Z2 vortices 
affects the temperature dependence of the ESR absorp- 
tion linewidth and the frequency dependence of the ESR 
resonance fields at 1.3 K below T* is well explained by a 
conventional spin wave theory. These results suggest an 
occurrence of Z<i vortex-induced topological transition. 

The ESR measurements for the frequencies up to 1.4 
THz were conducted in pulsed magnetic fields up to 
about 53 T at temperatures between 1.3 K and about 
90 K. The ESR measurements below 500 GHz in steady 
magnetic fields up to 14 T at 1.4 K were done by utiliz- 
ing a superconducting magnet and a vector network an- 
alyzer. Single crystals used for the measurements were 
grown by a chemical- vapour transport technique using 
iodine Q. 
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FIG. 1: Temperature dependence of ESR absorption spec- 
tra of NiGa 2 S 4 at 584.8 GHz for (a) H\\c and for (b) H±c. 
The arrow and cross indicate the broad large signal and the 
weak one, respectively. The sharp signal at about 21 T comes 
from an ESR standard sample of DPPH for correction of the 
magnetic field. 



Figures 1(a) and 1(b) show the temperature depen- 
dence of the ESR absorption spectra at 584.8 GHz in 
pulsed magnetic fields for i?||c and for if_Lc, respectively. 
We observed a strong resonance signal indicated by an 
arrow whose linewidth is extremely broad and the other 
weak one indicated by a cross. The weak resonance sig- 
nals are considered to come from the paramagnetic reso- 
nances caused by imperfection of the lattices. The tem- 
perature dependences of the absorption linewidth and the 
resonance field of the broad signal are plotted in Fig. 2 
and its inset, respectively. The broad ESR signal shifts 
toward lower field with decreasing temperature for both 
directions below about 30 K. The absorption linewidth 
increases as the temperature is lowered from about 80 
K to 8.5 K, and a bending appears at 23 K at which 
the increasing tendency of the linewidth with decreasing 
temperature becomes gradual. This behavior between 
80 K and 23 K is interpreted as the line broadening due 
to the development of the SRO as suggested from the 
NQR measurements Q. Thus, we separately examine 
the temperature ranges above 8.5 K, namely above 23 K 
and between 23 K and 8.5 K. First, the linewidth above 
23 K is investigated in terms of the critical behavior in 
the 2D Heisenberg AFM. The common expression of tem- 
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FIG. 2: Temperature dependence of full width at half maxi- 
mum of the ESR linewidth in NiGa 2 S 4 for H\\c at 584.8 GHz. 
Solid and broken lines indicate calculated line widths caused 
by the dynamics of Zi vortices and the development of 2D 
SRO, respectively. The inset shows the temperature depen- 
dence of the resonance field of the large ESR signal in NiGa2S4 
at 584.8 GHz. Open and closed circles represent the resonance 
fields for H\\c and HJlc, respectively. 



perature dependence of the linewidth (full width at half 
maximum) associated with the AF LRO temperature Tn 
is 



AH 1/2 <x(T-T N y 



(1) 



where p is a critical exponent, which reflects the 
anisotropy and the dimensionality of the system. Based 
on Kawasaki's dynamic equations, the critical exponent 
of 2D AFM is given by p — (3 — 277)1^, where v is the crit- 
ical exponent for the inverse correlation length and 77 the 
critical exponent related to the static correlations [ljj. 
In the case of 2D Heisenberg AFM, p is evaluated to be 
about 2.5 from the experimentally obtained r\ = 0.2 and 
v = 0.95 [16]. In NiGa2S4, one expects that the devel- 
opment of the SRO gives the same critical behavior as 
that of the 2D Heisenberg AFM. In addition, since there 
is no LRO down to 0.08 K Q, In must be small enough 
to approximate 0. Thus, we fit Eq. (1) with p = 2.5 and 
Tn — to the experimental line width below 80 K and 
obtained good agreement above 23 K. Below 23 K, the 
increase of the linewidth with decreasing temperature be- 
comes gradual. This indicates that some effect disturbs 
the development of the SRO below 23 K. We take into 
account a topologically stable point defect which is the 
so called Z 2 vortex, characterized by a two-valued topo- 
logical quantum number, suggested theoretically [17[ in 
the Heisenberg TLAFM. Unlike a conventional vortex of 
spins in the XY model, it is a vortex of spin chiralities. 
The Zi vortices form bound pairs at low temperatures 
and begin to dissociate at a certain critical temperature 
Ty, causing a Kosterlitz-Thouless-type phase transition. 
Above Ty the unbound free vortices are thermally ex- 
cited, and the spin fluctuates through the passage of the 



vortices. Hence, unbound Z^ vortices disturb the devel- 
opment of the SRO. The effect of Zi vortex on the ESR 
linewidth is discussed for Heisenberg TLAFMs, HCrC>2 
and LiCrC>2 [18| • According to the analysis used for these 
compounds, the ESR linewidth above TV is expected to 
behave as, 
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where Ey is an activation energy of a free Z2 vortex J18j 
and fee the Boltzmann constant. According to the es- 
timation by a Monte Carlo simulation [17|, Ty and Ey 
are given by T v = 0.31JS 2 and E v = 1.65JS 2 , where 
J is the nearest neighbor exchange interaction, and thus 
the ratio Ey /Ty becomes constant. We take TV = 8.5 K 
from the anomaly of the magnetic susceptibility, and then 
obtain Ey = 45.2 K. Using this value, satisfactory agree- 
ment between experiment and calculation is attained be- 
tween 8.5 K and 23 K as shown in Fig. 2. From the neu- 
tron scattering experiment, it was found that the correla- 
tion length in the triangular plane grows up to be about 
three lattice units at 23 K [lj|. This result indicates 
that the correlation length of about three lattice units 
is needed to observe the Zi vortex effect. The tempera- 
ture dependent shift of the resonance field below about 
30 K in the inset of Fig. 2 is also thought to be caused 
by growing up of the correlation length which is about 
two lattice units at 30 K [19(. Below about 10 K, the 
linewidth shows a gradual decrease as the temperature is 
lowered and even at the lowest temperature, it remains a 
large value. From the NQR measurements, Takeya et al. 
suggested that a gradual spin freezing occurs between 10 
K and 2 K, and inhomogeneous internal fields appears 
below 2K [|. This considerably large ESR linewidth at 
the lowest temperature must reflect the distribution of 
the internal fields. Accordingly, the temperature depen- 
dence of the ESR linewidth is consistent with the NQR 
results. 

Next, we discuss the frequency dependence of the res- 
onance field of the strong ESR signal at the lowest tem- 
perature. Figures 3(a) and 3(b) show the resonance fields 
at various frequencies obtained from the ESR absorption 
spectra at 1.3 K in pulsed magnetic fields and at 1.4 K in 
steady magnetic fields for H\\c and for H-Lc, respectively. 
The broken line in each figure represents a paramagnetic 
resonance line with g- factor measured at about 90 K. 
As mentioned before, the weak signal probably comes 
from isolated Ni ions due to the imperfection of the lat- 
tices and puts on the paramagnetic resonance line. The 
ESR modes of the strong signal have a zero-field gap 
and approach the paramagnetic resonance modes with 
increasing fields for both field directions. These ESR 
modes are reminiscent of spiral spin resonance mode ob- 
served in TLAFM with easy-plane anisotropy [2(|. It 
is known that even in the disordered phase, some low- 
dimensional magnets exhibit resonance modes expected 
in the LRO phase in the temperature region where the 
SRO well develops [2l|. Therefore, we analyze the exper- 
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FIG. 3: Frequency-field diagram of NiGa2S4 at 1.3 K in pulsed 
magnetic fields and at 1.4 K in steady magnetic fields for (a) 
H\\c and (b) H±c. Closed circles and crosses denote strong 
and week signals, respectively. Solid and broken lines repre- 
sent calculated resonance modes and paramagnetic resonance 
line, respectively. The insets in Figs. 3(a) and 3(b) show the 
calculated spin wave dispersion relation near k~0 at H=0 T 
where a is the lattice parameter and the magnetization curve 
for HJ-C at 1.3 K near the softening field around 44 T, re- 
spectively. 



imental results in terms of spiral spin resonances. A con- 
ventional mean-field approximation theory derives that 
the exchange constants are required to be J\j Jz=— 0.2 
and T2=0, where Ji<0 (ferromagnetic) and T3>0 (an- 
tiferromagnetic), to realize the 57° spiral spin structure 
indicated by the neutron diffraction measurements. The 
Ji (i — 1, 2, 3) represent the first, second, and third 
neighbor exchange interactions, respectively. The ex- 
istence of considerably strong third neighbor exchange 
interaction is also suggested in Ref. [T^]. Since J 3 is 
large enough compared to the others, we consider Ji=0 
for simplicity. From the conventional spin wave theory, 
three kinds of excitation modes are extracted for the spi- 
ral magnet. In the case of the Heisenberg system with 
no anisotropy, these modes are degenerate and gapless at 
zero field. However, when the system has the easy-plane 
anisotropy (T>>0), two modes have a gap whereas one 
mode remains gapless. These two excitations with a gap 
are expected to be observed in the ESR experiment. The 
resonance modes with 120° spiral spin structure for H\\c 



and iJ_Lc were derived by the conventional spin wave the- 
ory [22| and molecular field approximation [23|, respec- 
tively. The calculated three resonance modes, namely 
two gapped and one gapless modes, are drawn in Fig. 
3. The agreement between experiment and calculation is 
good by using the following parameter values; J3/fcB=21 
K, D/fc B =0.9 K, .g||=2.09, g±=2.15. The zero-field en- 
ergy gap is evaluated to be about 13 K. From the ob- 
tained J 3 , we estimate the Weiss temperature to be —84 
K, which is nearly equal to the value obtained from the 
magnetic susceptibility measurement [6| . We did not ob- 
serve the resonance modes, fux> 2 and fka^, which gradu- 
ally drop with increasing field. These modes have weak 
field dependence and are expected to show much broader 
resonance linewidth than the upper ones. Hence, these 
resonance modes were hardly observed. According to 
Ref. 23], the magnetic field at which the fkvs becomes 
zero corresponds to the occurrence field of the collincar 
spin configuration, in which spins on two sublattices are 
parallel and those on one sublattice are antiparallel to the 
applied magnetic field in the basal plane. Correspond- 
ingly, we observed an anomaly of the magnetization curve 
for H±c around 44 T as shown in the inset of Fig. 3(b), 
the magnitude of which is about one third of the satura- 
tion magnetization (2.15 /ie/Ni), while no anomaly was 
observed in the magnetization curve for H \\c. 

The above analysis suggests that spin-wave like exci- 
tation modes develop at low temperatures in NiGa2S4. 
Importantly, the lowest gapless mode derived from our 
analysis has fc-linear dispersion at fc ~ as shown in the 
inset of Fig. 3(a). This mode remains gapless even in fi- 
nite magnetic fields. Furthermore, spin wave velocity v s , 
which corresponds to a slope of the fc-linear dispersion, 
is almost constant in the field low enough compared with 
the saturation field, which is evaluated to be about 130 T. 
Actually, v s = 948 m/s at 7 T calculated for H \\ c with 
the parameters determined from the ESR measurements 
is almost the same as v s = 949 m/s at T. In addition, 
it is considered that the gapped modes does not affect 
the specific heat below 10 K because of the energy gap 
evaluated to be about 13 K. These features can qualita- 



tively explain the field independent specific heat with T 2 - 
dependence. From the magnon with fc-linear dispersion 
on the 2D triangular lattice, we obtain the specific heat at 
low temperature region as C — (3£(3) /ir)N \k^(T /v s ) 2 , 
where Na is the Avogadro number and ((3) — 1.202. 
The calculated result, however, is somewhat smaller than 
the experimental one. For a quantitative agreement with 
the experimental specific heat, the value of v s is required 
to be nearly three times smaller than that expected from 
our ESR analysis. From the conventional spin wave the- 
ory, a well-defined sharp energy dispersion is expected. 
On the other hand, the distribution of the internal field in 
NiGa2S4 likely gives broad energy structure of the exci- 
tation modes, which were implied from the considerably 
large ESR linewidth. This broadening of the energy spec- 
trum may cause specific heat larger than that expected 
from the spin wave theory. For direct observation of the 
energy dispersion, detailed inelastic neutron scattering 
experiment is desired. 

Finally, we discuss the possibility of the Z 2 vortex sce- 
nario in Ref. [XOl ] . If the Z 2 vortex-induced topological 
transition with small biquadratic exchange interaction 
occurs at T v , the low-temperature phase is dominated 
by spin wave excitations, and the spin correlation length 
remains finite [10]. In our observation, below T v ~8.5 K, 
the experimental results suggest spin wave like excita- 
tions. In addition, the Z 2 vortices effect can be found 
in the temperature dependence of the ESR absorption 
linewidth. Consequently, the unconventional properties 
of NiGa 2 S4 in the low temperature region strongly sug- 
gests such a vortex-induced topological transition. 
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